Boltzmann theory of magnetoresistance due to a spin spiral 
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We studied the magnetoresistance due to a spin spiral by solving the Boltzmann equation. The 
scattering rates of conduction electrons are calculated by using the non-perturbative wave function 
of the conduction electrons and the non-equilibrium distribution function is obtained by numerically 
solving the Boltzmann equation. These enable us to calculate the resistivity of a sufficiently thin 
spin spiral. A magnetoresistance ratio of more than 50 % is predicted for a spin spiral with high 
spin polarization (>0.8) and a small period (about 1-2 nm). 
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There is great interest currently in spin-dependent 
transport phenomena in magnetic domain walls such 
as the magnetoresistance (MR) efFectii^i^iiii^ and spin- 
transfer torque-driven magnetization dynamics^'i'^'^ be- 
cause of the potential application of these phenom- 
ena to spin-electronics devices such as spin-motive-force 
memorjiiS and racetrack memory^^. In these devices, 
higher magnetoresistance due to a thin domain wall is 
desirable for high-density magnetic recording. 

In 1997, Levy and Zhang^ studied the resistivity due 
to domain wall scattering by using the same Hamiltonian 
that was used to explain the giant magnetoresistance ef- 
fect. They found that the magnetoresistance ratio is pro- 
portional to where d is the thickness of the domain 
wall, and showed that the magnetoresistance ratio is be- 
tween 2% and 11%, which is consistent with the exper- 
imental results (5%) of Refi where the thickness of the 
domain wall is about 15 nm. 

However, the theory of Levy and Zhang^ cannot be 
applied to a sufficiently thin domain wall for two rea- 
sons. First, the scattering rates of the conduction elec- 
trons are calculated by using the pcrturbative wave func- 
tion, which is up to the first order of the dimensionless 
parameter ^. The parameter ^ = Ij/d characterizes the 
non-adiabaticity of the spins of the conduction electrons 
with respect to the localized spins, where Ij = i:hv-p / {'^J) 
is the electrons' traveling length during the precession of 
their spins around the sd-exchange field J. For a do- 
main wall with ^ > 1, the theory cannot estimate the 
amount of the non-adiabaticity correctly, and thus can- 
not be applied. Second, since Levy and Zhang applied 
the diffusion approximation to the Boltzmann equation, 
their theory cannot be applied to the domain wall in the 
ballistic region d<lraip^ where Imip is the mean free path. 
For conventional ferromagnetic metals such as Fe, Co, 
Ni, and their alloys, both Ij and ^mfp are on the order of 
a few nm^^. 

The thickness of a domain wall is determined by the 
competition of the exchange coupling between the local- 
ized magnetizations and the magnetic anisotropy, and 
is usually on the order of 50 nm for conventional ferro- 
magnetic metals. Recently, however, the production of 



the domain wall of Co5oFe5o, with a thickness of about 
2.5 nm, was achieved by trapping the domain wall in a 
current-confined-path (CCP) geometry^^, and a magne- 
toresistance ratio of about 7% -10% was observed. Many 
studies have examined to understand the physical prop- 
erties of the CCP structure and applied that structure to 
magnetic device&i^i-^^ . To investigate the transport prop- 
erties of such a thin magnetic structure, in which the 
system size d is comparable to or less than Ij and Zmfp, 
i.e., a few nm, it is important to develop the theory of 
Levy and Zhang to take into account the amount of the 
non-adiabaticity correctly and to describe the transport 
without the diffusion approximation. 

In this paper, we study the dependence of the magne- 
toresistance ratio of a spin spiral on its period (thickness) 
d by solving the Boltzmann equation. We extend the the- 
ory of Levy and Zhang^ by using the non-perturbative 
wave function of the conduction electrons in the calcula- 
tion of the scattering rates and by solving the Boltzmann 
equation of the non-equilibrium distribution function nu- 
merically. These enable us to investigate the resistivity 
due to a spin spiral with d<lj, Imip- We find that the MR 
ratio is more than 50% for a spin spiral with high spin 
polarization (/?>0.8) and a small period (d ~ 1 — 2nm). 
We also find that in the diffusive region, d>lj ,lraip, the 
MR ratio is proportional to 1/c?^, while in the ballistic re- 
gion, d<?,/,^mfp, the MR ratio increases with decreasing 
d more slowly than it does in the diffusive region. 

We consider electron transport in a one-dimensional 
spin spiral that lies over —d/2 < x < d/2, where d is 
the period of the 7r-rotation of the localized spins. We 
assume that the spin-dependent transport of the conduc- 
tion electrons is described by the following Hamiltonian: 



Ho = - — V^-ja-S{r), 
Zm 



(1) 



where J is the sc?-exchange coupling constant between 
the conduction (s-like) electrons and localized (d-like) 
spin, (T is the vector of the Pauli matrices and S = 
(0, — sin 61, cos 6*) is the unit vector along the direction 
of the localized spin. The angle 9 is given by d{x) — 
{'K/d){x-\-d/2). On the other hand, the spin-dependent 
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impurity scattering is described by^ 

^^E^-.^'^-SWj'^lr-RO, (2) 

i 

where Ri is the position of the impurity, and v and j 
are the spin-independent and spin-dependent scattering 
potentials, respectively. The dependence of the transport 
properties on the direction of the electrons' spin arises 
from either the exchange energy J or the spin-dependent 
scattering potential j, i.e., the spin dependence of the 
number of the conduction electrons at Fermi level is due 
to J, and the spin dependence of the scattering rate is 
due to j. 

The resistivity of the spin spiral is calculated by solving 
the Boltzmann equation of the non-equilibrium distribu- 
tion function /^(k) given by 



-ewJi?(5(eF-e(k, s)) 



d^k' 
(2^ 
d^k' 



w^kk'[r(k)-r(k')] 



(27r) 



[r(k)-/-(k')], 

(3) 



where W^l^, is the scattering rate of the conduction elec- 
trons from the state (k, s) to the state (k', s'), £f is the 
Fermi energy and E is the strength of the applied electric 
field. The index s, s' = ± denotes the eigenstate of Hq in 
spin space, which is given byi^ 
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Here the angle (f>{kx) and the spinor rjj. are given by 
<f>{kx 



- arctan 



(5) 



V- 



(6) 



where 9' = d9/dx = TT/d and kj = \/2mJ/h, respectively. 
The factor tan((/)/2) characterizes the non-adiabaticity 
of the spins of the conduction electrons with respect 
to the localized spins, and is the most important pa- 
rameter in our calculations. It should be noted that 
this factor is always less than unity for any period d 
and momentum kx. For a sufficiently large period d, 
tan((/)/2) {kx0')/{2kj) = (kx/kF)^, and the wave func- 
tion ^ is reduced to the wave function calculated by 
Levy and Zhang**. On the other hand, for a small period 
d where ^ — Ij/d is comparable to or larger than miity, 
the wave function ^ does not equal the wave function 
given in Ref4.. The eigenvalue of Hq is given by 



£(k,s) 



2m 



'.^{kx9'y+k^j 



(7) 



The velocity is given by v^ = de{k, s)/dpx- The scat- 
tering rates are calculated by using the Fermi golden rule 
with the Born approximation. 
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W^kk' = Yl^kk;i''5(£(k,s)-e(k',s')), 



(8) 



where the matrix elements of the sattering potential ^ 
are calculated by using the wave function ^ and are 
given by 



IK 



kk' 



(v — sj) COS — cos \-(v + si) sin — sin — 

\ jj 2 2 2 2 
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{—sv+j) COS — sill — + {sv+j) sin — cos — 



(10) 

respectively, where Ci is the impurity concentration. 
Here, for simplicity, we denote (f>{kx) and 4>{k'^) as (j) 
and (/)', respectively. In the limit of d cxd, the con- 
duction electrons change the direction of their spins adi- 
abatically, and thus, tan((/)/2) — > for any momentum 
kx- In this limit, the spin-flip scattering rate is zero, 
i.e., V-^^,'^ = 0, and the spin-conserved scattering rate, 
^kk' °^ l^kk'P' independent of the momentum kx- 
On the other hand, in the limit of d— > 0, tan(0/2) 1 
for the large momentum kx — kp, which means that the 
amount of non-adiabaticity is maximized for the conduc- 
tion electrons with ~ vf because the traveling time 
through the spin spiral of these electrons, d/vx, is much 
shorter than the period of the precession of the spins of 
the conduction electrons around the exchange field J. In 
Ref,- , Levy and Zhang approximate that cos((/)/2) 1 
and sin(0/2) ^ tan(0/2) ^ (fcaj/Zcp)?- It should be noted 
that for a thin spin spiral where S, = lj/d is comparable 
to or larger than unity, the estimation of the scattering 
rate W^^, in our theory for large momentum kx is much 
smaller than that obtained by Levy and Zhang because 
the factor tan(0/2) in our calculation is always less than 
unity while the factor {kx/kF)S, used in Refi^ is larger 
than unity. Since the resistivity is high for a high scat- 
tering rate, the magnetoresistance obtained in our theory 
is much lower than that obtained by Levy and Zhang, as 
shown below. 

To obtain the non-equilibrium distribution function 
/'' (k) from the Boltzmann equation ([3]), we assume that 
r(k) = (9/^(0) (k)/ae)g^(k) ^ ~6{eF ~ e{k,s))g^{k), 
where /"(o) (k) is the distribution function in equilibrium. 
Then, Eq. ([3]) is reduced to 
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wave vector [nm"^] wave vedor [nm"^] 

FIG. 1: The dependence of the distribution function, /eE, 
on the momentum for (a) d — Inm and (b) d — lOnm, 
respectively. 

where kp is given by 



fcF = yfc|+(|) +s^{kF9')^ + kj . (12) 

The relaxation time T^{kx) is given by \/T^{kx) = 
l/r^^lkx) + 1/T^~^{kx), where the spin-conserved re- 
laxation time T^'^{kx) and the spin-flip relaxation time 
T'^~^{kx) are given by 

^^^^T'^^^l^kk'?- (13) 

The distribution function /''(k) is obtained by numer- 
ically solving Eq. (fTTj) ^^. The resistivity of the spin 
spiral is calculated as p = l/(fT++(T^), where cr* = 
-{e/E) J d^k/{2Trfv'^f'{'k) is the conductivity of the 
spin-s electrons. 

In the calculation of the scattering-in term, 

/d3k'/(2^)3[iy-,/«(k') + Ty^k'7-^(k')], in Eq. m, 
Levy and Zhang"* assume that the non-equilibrium dis- 
tribution function is proportional to the momentum k^- 
However, we do not apply this diffusion approximation 
to the scattering-in term because we are interested in the 
resistivity for a spin spiral with d< Imfp- Figure[T](a) and 
(b) show typical dependences of the distribution function 
obtained by Eq. /eE, on thc momentum kx 

for d = 1 nm and = 10 nm, respectively, where the 
mean free path l-adp is taken to be 5.9 nm. According to 
Fig. [Jl we can verify that the diffusion approximation is 
not applicable to the region d < l^fp while it is a good 
approximation to the region d>l^[p. 

Before estimating the resistivity of a spin spiral, we 
should emphasize the validity of our calculation. The 
semi-classical Boltzmann equation is applicable when the 
system is larger than the width of the wave packet of 
the conduction electrons, i.e., the Fermi wavelength Ap. 
In our calculation, this condition equals c? > Ap. For 
conventional ferromagnetic metals, the Fermi wavelength 
is on thc order of a few angstrom, which is one order of 
magnitude smaller than Ij and Imf^- It should also be 
noted that the derivative of the angle 9{x) is assumed to 
be constant in the derivation of the wave function Q. 
Thus, our calculation is valid for a spin spiral where the 
direction of the localized spin changes linearly in space. 



Figure [2] shows the dependence of the MR ratio due 
to a spin spiral, defined by {p — p^^'')/ p^^\ on its period 
d. The values of the parameters we use are as follows. 
The Fermi energy ep and the sd-exchange coupling con- 
stant J are taken to be 5.0 eV and 0.5 eV, respectively. 
The Fermi wavelength Ap is estimated to be 5.4 A. The 
strengths of the impurity scattering, v and j, and the 
impurity concentration, q, are estimated by the resis- 
tivity p*-*^) and the spin polarization /3 of a bulk ferro- 
magnetic metal. The value of is taken to be 150 
rinm, which is a typical value of the conventional ferro- 
magnetic metals while the value of /3 is taken to be 
from 0.3 to 0.9. Using these parameters, Ij = ttHvf/{4:J) 
is estimated to be 1.4 nm, and the mean free path 

Up = (Cfp+^mfp)/2, where /^fp = 4^^'"\ ^| = fifcfV m, 
T«(o) =7r?iV[mci(t;-sj)2fc;(°)], and fc;(°^ = ^k'^ + sk'j, is 
estimated to be 5.9 nm, which is approximately indepen- 
dent of the values of (3. 

As shown in Fig. [21 the MR ratio increases as the pe- 
riod d decreases. The higher the spin polarization of the 
bulk /3 is, thc higher the MR ratio is. In the diffusive re- 
gion d> /j, ^mfp, the MR ratio is estimated to be l%-20%. 
On the other hand, for a thin spin spiral (d ~ 1 — 2nm) 
with a high polarization (/? ~ 0.8 — 0.9), an MR ratio 
of more than 50% is predicted. Recently, a spin spiral 
of ferromagnetic Mn/W(001) with the rotation period 
2d ~ 2.2 nm was created experimentallj^, whose period 
d is comparable to or smaller than Ij and /mfp- Thus, it 
is reasonable to consider such a sufficiently thin spin spi- 
ral d<lj,lynfp- The values of the spin polarization /3 of 
the conventional ferromagnetic metals such as Fe, Co, Ni, 
and their alloys are about 0.5-0.7; for example, /3 = 0.51 
for Co, 0.65 for CogiFcg, and 0.73 for NigoFeac,^'^. The 
value of /3 depends on the combination and the composi- 
tion ratio of the ferromagnetic metals, and we can expect 
ferromagnetic metals with high spin polarizations. Thus, 
the prediction of our calculation for a spin spiral with 
high spin polarization /3 and a small period d < Ij, Zmfp 
will be confirmed experimentally. 

The physics behind these results are as follows. The 
origin of MR due to a spin spiral is the mixing of the chan- 
nels of the spin-up current and spin-down current due to 
the spin-dependent scattering potential V. The channel 
mixing increases the scattering probability of the conduc- 
tion electrons, and thus the resistivity. The mixing due 
to the scattering arises from the non-adiabaticity of the 
spins of the conduction electrons, which is characterized 
by ta,n[(l){kx)/2]. In the limit of d^oo, the conduction 
electrons change the direction of their spins adiabatically, 
i.e., tan((/)/2)— s-0 for any momentum k^, and the MR ra- 
tio tends to be zero. On the other hand, in the limit of 
d—^O, the amount of non-adiabaticity that is maximized 
for the conduction electrons with large momentum k^, 
i.e., tan((/)/2) 1 for k^ — kp, and thus the MR ratio, 
increase as the period d decreases. In other words, the 
MR due to the spin spiral is mainly due to the conduc- 
tion electrons with large momentum kx- Since the MR 
arises from the asymmetry of the transport properties of 
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period of spin spirai, d [nm] 

FIG. 2: The dependence of the magnetoresistance (MR) ratio 
of a spin spiral on its period d. The sohd hnes from bottom 
to top correspond to the MR ratio with the spin polarizations 
/?=0.3, 0.4, 0.5, 0.6, 0.7, 0.8 and 0.9, respectively. The dashed 
line is the MR ratio estimated by the theory of Levy and 
Zhang* with /3 = 0.5. 

the spin channels, the higher the spin polarization /3 is, 
the higher the MR ratio is. 

The dashed line in Fig. [2]shows the MR ratio estimated 
by the theory of Levy and Zhang with /3 = 0.5^; 

By comparing the solid line and the dashed line in Fig. 
[21 we find that the MR ratio in the diffusive region, 
d > Ijjlmfp, is proportional to as shown by Levy 

and Zhangi. On the other hand, in the ballistic region, 
d < Ijjlmip, the MR ratio increases more slowly as the 
period d decreases compared to the diffusive region. It 
should be noted that the factor tan(0/2) is approximated 



to be (kx/kp)^ in Ref.^, which is on the first order of 
However, for a thin spin spiral, the higher-order terms of 
1 /d also contribute to the calculations of resistivity, and 
the dependence of the MR ratio on the period d shifts 
from As shown in Fig. [2l the MR ratio obtained 

by our theory is much smaller than that obtained by 
Levy and Zhang. This is due to the fact that the esti- 
mated scattering rate by our calculation is much lower 
than that by Levy and Zhang, as mentioned above. The 
smaller the period d is, the larger the difference is in the 
amount of non-adiabaticity between our theory and that 
of Levy and Zhang, i.e., the difference in the values of 
tan((/)/2) and {kx/k-p)£^. Thus, the difference in the MR 
ratio between our theory and theirs increases as the pe- 
riod d decreases. 

In conclusion, we have studied the dependence of mag- 
netoresistance due to a spin spiral on its period d by solv- 
ing the Boltzmann equation. The scattering rate of the 
conduction electrons in the spin spiral is calculated by us- 
ing the non-perturbative wave function of the conduction 
electrons, and the non-equilibrium distribution function 
is obtained by numerically solving the Boltzmann equa- 
tion. An MR ratio of more than 50% is predicted for a 
thin spin spiral (d~ 1 — 2nm) with high spin polarization 
(/3>0.8). We also find that the MR ratio in the diffusive 
region is proportional to l/d^, while in the ballistic re- 
gion the MR ratio increases more slowly with decreasing 
d compared to the diffusive region. 
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